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Abstract. We consider small perturbations of the Laplace operator in a 
multi-dimensional cylindrical domain by second order differential operators 
with periodic coefficients. We show that under certain non-degeneracy condi- 
tions such perturbations can open a gap in the continuous spectrum and give 
the leading asymptotic terms for the gap edges. We also estimate the values of 
quasi-momentum at which the spectrum edges are attained. The general ma- 
chinery is illustrated by several new examples in two- and three-dimensional 
structures. 



1. Introduction 



& ■ This paper is devoted to the spectral problem for a class of periodic differential 

operators in unbounded periodic domains. The spectra of such operators are known 
to have band structure, and one is interested in the question whether the spectrum 
fills a semi-axis or has gaps. The importance of such questions is motivated, in 
particular, by the study of photonic crystals [JJ: the spectral gaps correspond to 
the energies at which no photons can be transmitted through the material. 

Existence of gaps for various periodic operators attracted a lot of attention in 
the last decade, see e.g. the reviews We mention, in particular, the recent 

papers [4j [5j [6] discussing gap opening for tube-type domains with various pertur- 
bations. Our interest to the problem was attracted by the paper [7J [5], in which 
the authors discussed the locations of the extrema of the band functions inside 
the Brillouin zone for generic periodic operators. In particular, the paper [5] was 
devoted specifically to the operators with the one-dimensional periodicity, and it 
was shown that the spectral edges are generically different from the periodic and 
the anti-periodic eigenvalues as it is sometimes assumed. One of the most natural 
classes of one-periodic systems is provided by the waveguides, but it seems that 
there are just few works addressing the related questions. The paper [9J discussed 
the spectra of acoustic structures, and our papers [HHHI] showed the occurrence of 
the same effect for a configuration consisting of two interacting waveguides modeled 
by the free Laplacian, and this discussion was continued in the very recent work [12j 
at the example of a cylinder with a periodic system of holes. We mention another 
recent paper [TJ5] discussing related questions at the physical level of rigor. 

In the present paper we prove a general result giving sufficient conditions for the 
gap opening for a class of second-order operators in multi-dimensional cylindrical 
domains. It is shown that the presence of gaps at the energies different from the 
periodic and anti-periodic eigenvalues is a generic fact, and we discuss the parame- 
ters controlling the gap opening at various values of the quasi-momentum. Like in 
[Till ITT1 [T2] , the gap opening may occur at the points where the band functions of 
the reference operator (the free Laplacian) meet in some specific way, but the fur- 
ther analysis is done with the help of the elementary tools of the regular first-order 
perturbation theory, while the previous works used singular asymptotic expansions 
or other methods which are rather sensible to the type of perturbation. We discuss 
in detail several new examples and show that the gaps at the interior points of 
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the Brillouin zone and at its center are controlled by different parameters. In par- 
ticular, we put in evidence several perturbations which open gaps at intermediate 
values of the quasimomentum, but the leading asymptotic term appears to vanish 
at the center of the Brillouin zone. Therefore, our analysis considerably extends 
the family of situations for which one can prove the existence of gaps compared to 
the previous works. 

2. Formulation of the problem and the main result 

Let x' — (xi, . . . , x n -i), x — (x 1 , x n ) be Cartesian coordinates in R" -1 and R n , 
respectively n ^ 2, and u be a connected bounded domain in R™ _1 with Lips- 
chitz boundary. By ft we denote an infinite straight cylinder in R" with the base 
uj, namely, J] := w x R, and by Ho we denote the positive Dirichlet Laplacian in 
1/2 (fi) on the domain W% (tyt which is introduced as the subspace of the functions 
in W%(Q) with zero trace on dft. We assume that the boundary dco is sufficiently 
regular for the operator Ho to be self-adjoint in Lz^l). Introduce another differen- 
tial operator C e with the domain W$(£l) whose coefficients may depend on a small 
positive parameter e, 

i.j — 1 j — 1 

where Aij € W^(H.) are complex- valued functions satisfying Aij = Aji, and Aj € 
W^fi), Ao € Loo(ri) are real- valued functions. All the coefficients are supposed 
to be T-periodic in x n : 

Aij {x , x n + T, e) = Aij {x, e), i,j = l,...,n, 

Aj(x',x n + T,e) = Aj(x,s), j = 0,...,n, 

with an e-independent period T. We assume that for e — > +0 one has 

n n 

^\\A ij (.,e)-A ij (-,0)\\ wio{n) + ^2\\A j (.,s)-A j (-,0)\\ wUn) 
i,j=l °° j=i (2.1) 

+ \\A a (-,e) - A (-,0)|| Loo(n) =: 77(e) -> +0, 

where □ := u> X (0,T) is the elementary cell. The condition (|2.ip is satisfied, if, for 
example, all the coefficients are continuously differentiable in both x and e. 

Under the above assumptions and for sufficiently small e the operator C e is 
symmetric and relatively bounded w.r.t. Ho, and, by Kato-Rellich theorem, for 
small e one can define the sum H £ := Ho + eC e which is a self-adjoint operator 
in I/2 (^)- One can consider H e as the Hamiltonian of a quantum particle in the 
waveguide ft. 

Due to the periodicity, the operator H e can be studied using the Floquet-Bloch 
theory [TJ. Introduce self-adjoint operators 

H e (r) :=-A x ,+ +e£ e (r), (2.2) 

n n 

£ e (r) := £ 2i(r)^-/j(T) + £ (A^r) + l^A,-) + A , (2.3) 
jj=i 3=1 

/ ^ T ' ):=i ^' ^' = 1 '---' n_1 ' /n ^ :=i ^ _T ' 

depending on the parameter r € S = (— tt/T, tt/T] and acting in L2(0) on the 
domain W% per {0) consisting of the functions from W% (□) with zero trace on dco x 
(0, T) and satisfying periodic boundary condition at x n = and x n = T. The 
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parameter t is referred to as the quasi-momentum and the set B is usually called 
the Brillouin zone. 

For any r, the operator H £ (t) has a compact resolvent, and its spectrum consists 
of an infinite sequence of discrete real eigenvalues E m {r,£), m ^ 1, which we 
assume to be ordered in the ascending order counting multiplicity. The mappings 
t i— > E m (r,s) are called band functions; it is known that they are continuous with 
E m (~Tr/T,e) = E rn (Tt/T, e) and that the spectrum of H E is the union of their 
ranges, 



oo 

<rCH E )= (J {E m (T,e) : r e s}. 



In particular, mia{'H £ ) = S e = min re g E\(t, e). Each open interval (a,/3) C 
[£ £ ,+oo) \ o-(W e ) with a, P £ o-(W e ) is called a gap of "H e . We note that the 
unperturbed operator T-Lq has no gaps and its spectrum fills a semi-axis (see below). 
We also observe that the band functions are not constant on any interval, and the 
spectrum of H e is absolutely continuous HH1 IIZ1 HB] ■ The main aim of the present 
paper is to obtain some conditions guaranteeing the existence of gaps for T-L e and to 
determine the values of the quasi-momentum r at which the band functions attain 
the respective gap edges a and j3. 

The eigenvalues of Hq(t) can be found by the separation of variables. Denote 
by —A.u the positie Dirichlet Laplacian in u). Since this operator is self-adjoint 
and has a compact resolvent, its spectrum consists of real eigenvalues of finite 
multiplicity denoted by Xj, j — 1,2,..., and assumed to be ordered in the ascending 
order counting multiplicity. The spectrum of Ho (r) then consists of the eigenvalues 

Ai J ,(r)=A i +^r+^ , j e N, P eZ, (2.4) 

and the values E m (r, 0) are obtained by their rearrangement in the ascending order. 
Moreover, if ipj(x') are the eigenfunctions of — AL associated with the eigenvalues 
Xj and chosen real and orthonormal, then the eigenfunctions ^ j :P of Hq(t) for the 
eigenvalues Aj -P (r) can be written as 

^W^^*^ 1 ". (2-5) 
and they are orthonormal in Let us formulate our main result. 



Theorem 2.1. Let numbers To € 
following conditions: 



7T ' 

°'r. 



j, k € {1,2} and p,q € Z satisfy the 

Xj and Xk are simple eigenvalues of— > (2-6) 

A i)P (r ) = A fe , g (r ) =: A ; (2.7) 

(ro)^(roXO; (2.8) 

A ^ A z , s ([0,^]) as(l,s)t{(j,p),(k,q)}; (2.9) 

(jC (±To)^j, p , *fc,,) La (n) ^ 0. (2.10) 

Define the functions 

P±(r) := ± J|^l^(r)- fc ?(r)- + fc 2 (r), (2.11) 
F3(r)| v fe 3 (r) 
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where 



(2.12) 



fc 3 (r) := -{p-q), fc 4 (r) : 

So(r) 



T r v 7 2 

b°n(r) &$ 2 (r)\ /X4>to*Ap. (A^t)^,*^^ 
^iW b° 22 (r)J l v (£o(r)* fe , 9 ,^- p ) i2(n) (£ (r)* fc , g ,* fc , 9 > i2(n)y 



(2.13) 



(2.14) 



(2.16) 



t ^ 0, and 

njp + q) bUr)+bl 2 {T) 
k i{T) := r, fc 2 (r) := , 

fe(r) := -(g-p), fc 4 (r) := - 11 w . 

So(T) -UiW tSaWj (2 ' 15) 
= /<A(r)^-,-p,^,-p> i2(n) <A(r)* J -,-p,*fc ) - 9 > £a(P A 
^(A)(r)* fcl -„*i,-p> £a(n) (£ (T)$ l ,- gi $ trt ) ia(D) j 

as t < 0. TTien 

(Al) XTie strict inequalities /3_(to) < /3 + (ro) arcrf T o) < /?+( — r o) ftoZd true. 
(A2) //, in addition, one has 

A := max{/3_(r ),/3_(-r )} < /3 r := min{/3 + (r ), /3+(-r )}, (2.17) 

iften ttere exists Eq > swc/i i/iai /or aZ/ £ € (0, £q) fie operator H, e has a 
spectral gap (ai(e), a r (e)) whose edges have the asymptotics 

a l/r (e) = A +e(3 l/r + O(e 2 +e V (e)), (2.18) 

anrf i/ie associated band functions Ei/ r {e,r) of H £ attain the gap edges 
ai/ r (e) at the points Ti/ r (e), 

mmE r (e,T) = E r (e,T r (e)), win Ei{e,t) = ^(e, 77(e)), (2.19) 

r r 

/or which the asymptotics 

n/ r (e) = T + £ 7i/r + 0(e 3/2 + en 112 ) (2.20) 
fc i( T *)l & i2( T *)l Mr*) 



7i 



fei(r.)|62 2 (r,)| h(r*] 



|fc 3 (n)| V^(r*) 2 - A?(r„) ks(r,)' 

are valid. In each of the latter formulas the number r* should be chosen to 
that of the values ±tq, ai which the maximum or minimum is attained in 
the formulas {E.lty . 
(A3) If at least one the following two conditions is valid: 

• all the coefficients H e are real 

• t = 0, 

then the condition is satisfied. 

Before proceeding to the proof, let us give some comments on the assumptions 
and the statement of the theorem. 
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Remark 2.1. Eq. (|2.6[) is a non-degeneracy condition. It is independent of the 
period T and concerns only the eigenvalues of the cross-section operator — A« . 
This condition allows us to reduce the spectral study of H e in a vicinity of Ao to an 
eigenvalue splitting problem for a 2 x 2 matrix. Note that this condition is always 
valid for j = k = 1, because the lowest eigenvalue of the Dirichlet Laplacian in a 
bounded connected domain is always non-degenerate. 

Remark 2.2. The conditions (|2.7|) . (|2.8[) and (|2.9I) impose some restrictions on 
the behavior of the band functions Aj )P and A& jg an the intersection point. The 
condition (|2.7p means exactly the presence of an intersection. The condition (I2.8[) 
shows that one of the functions increases and the other decreases at the intersection 
point. Finally, the third condition (|2.9p expresses the fact that the projection of the 
intersection point on the ordinate axis should not be overlapped by the projections 
of the graphs of the remaining band functions A;. s (r) as r g [0,n/T]. We observe 
that it is possible only as j + k ^ 3. 

One can easily see that the validity of these assumptions is conditioned by the 
presence of certain relations between the transversal eigenvalues Aj and the period 
T. We remark first that the conditions (|2.4I) and (|2.7I) can hold true only for p and 
q having opposite signs. Moreover, due to (|2.4[) one has 

A fc -A, = ^^(^±^ + 2.„). (2.21) 

In particular, if all three conditions in question can be satisfied with tq = for 
j = k = 1 and q = —p only, and if all the conditions are valid for To = tt/T, then 
automatically j = k = 1 and q = —p — 1. 

To guarantee the presence of at least one combination (fc, q) satisfying the 

above assumptions for To = 0, it is sufficient to ask for the inequality Ai i i(0) < 
A2,o(0), which is equivalent to T > 2ir/ — Ai; in this case all three conditions 
are valid for To = 0, j = k = 1, p = 1, q = —1. On the other hand, to satisfy 
the conditions (|2.7[) and (|2.8p for some t £ (0,n/T) it is sufficient to obey, for 
instance, the inequality A2.o(0) < Ai i i(0), then one can take j = 1, k = 2, p = — 1 
and q = 0, and the condition (|2.9[) is equivalent to A2,o('?o) < A3 i o(0). Rewriting 
these inequalities with the help of the explicit expressions (|2.4[) . one can easily see 
that the three conditions (|2T7)l . (|2~%|t . (pH)]) hold true with some t <E (0, n/T), j = 1, 
k — 2, p — — 1 and q = 0, if the period T obeys the estimates 

2tt 2tt 
VA3 — Ai — VA3 — A2 VA2 — Ai 

With the help of the explicit formulas for Aj tP one can easily construct other suf- 
ficient conditions. Note that for any fixed cross-section uj one can satisfy the as- 
sumption by an appropriate choice of the period T, and, moreover, one can obtain 
in this way any prescribed finite number of intersections satisfying all the assump- 
tions. One is also able to choose parameters in such a way that the assumption will 
be satisfied for any prescribed value of To- 

Remark 2.3. The key condition (I2.10j) allows one to make some conclusions on the 
behavior of the perturbed band functions in a vicinity of To using some elementary 
tools of the regular perturbation theory. Checking this condition is usually the 
most difficult step when applying the theorem, which will be seen with the series of 
examples below. As follows from IA3l for operators with real- valued coefficients as 
well as for To = 0, Eq. (12.10)) is the only non-trivial assumption on the perturbation 
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Remark 2.4. The theorem can be extended to other types of unperturbed operators 
as far as the unperturbed band functions have a similar structure. We do not de- 
velop this direction to avoid technicalities. We just remark that the results are valid 
in literally the same form for some other boundary conditions dQ, in particular, for 
the Neumann condition and the Robin condition with a constant coefficient, and 
such a modification results in an appropriate redefinition of the eigenvalues Xj and 
the eigenfunctions tpj. 

3. Proof of the main result 

This section is completely devoted to the proof of Theorem 12.11 We do all the 
constructions for r £ [0, ir/T) only; the study of r = n/T is completely identical, 
but it requires cumbersome notation since it corresponds to the end point of the 
Brillouin zone. 

By (|2.1[) one can represent C e = Co + r](s)A4 e , where M. e the operator is of 
the same type as C £ except the fact that instead of the continuity in e at e = 
we ask for the uniform boundedness for e < E\ in the appropriate Sobolev norms. 
Under these assumptions and for sufficiently small e the operator M. e is relatively 
bounded w.r.t. JC e :— Hq + eCq. By analogy with (|2.3p one can define the operators 
M e {t) and /C e (r), and then M e (t) is also relatively bounded w.r.t. fC e (r). 

Denote by E m (r, e), m = 1,2,..., the eigenvalues of JC e (t) taken in the ascending 
order counting multiplicity. By the minimax principle one can choose a constant 
c > such that for sufficiently small e the uniform in r £ B estimates 

\E m ( T ,E)-E m (T,0)\ <C£|£ m (r,0)|, \E m (T,E)-E m ( T ,e)\ ^ cerf{e)\E m {r, 0)| 

(3-1) 

hold true. Since the required estimates in lA2l are linear in e, the second inequality in 
(|3.ip shows that it is sufficient to consider the case of an e-independent perturbation, 
C e = Co, which will be assumed throughout the rest of the proof. 

The explicit expressions show that for any C > Ao there exists just a finite 
number of the pairs (I,s)eNxZ with (-co, C) n A; )S (B) ^ 0. Pick an arbitrary 
C\ > Ao, then, by (|3.1[) . one can find N £ N and £2 > such that E m {r,E) > 
C\ for all m > N, e < £2 and t £ B. Thus, there exists C2 > such that 
\E m (r, e) - E m (r, 0)| < C 2 e for all m < N, \e\ < e 2 and t £ B. 

Now we find M £ N for which A = Em(tq, 0) = Em+\{tq, 0). By the preceding 
estimates we can conclude that for any 8 > there exists £3 > such that (Ao — 
6, A + 5) n Ej {B, e) = for all j £ {M, M + 1} |e| < e 3 . Thus, the spectrum of U e 
near Ao for small £ is determined by the behavior of the two band functions Em 
and Em+i near ±ro, and for e — they coincide locally with the functions Aj tP 
and A k>q . 

By (|2.4p we have Aj i _ p (— r) = A J:P (r). This shows that the assumptions (12. T[) . 
(|2.8|) and (|2. 10|) are also valid with j, p, k, q and To replaced by j, —p, k, —q and 
—To, respectively, and with the same value Ao. One can find a neighborhood of 
Ao having no intersections with the ranges of A;, s as (l,s) £ {(j,p),(k,q)} and, 
moreover, it is clear that in [— n/T, tv/T] there are no other values ±to for which 
Eq. (|2.7p holds. To summarize, for £ £ (0, £3) and for any a > there is a constant 
to(a) > such that 

dist (a (H e (T)), A ) > ae for \t ± r | > t Q (a)e. (3.2) 

In other words, for the indicated values of r the spectrum oi% e (r) is separated from 
Ao by a distance at least Ae, and it is now sufficient to study the band functions 
of 'He near A as |r ± tq\ < t (a)e. We consider the case of To only; the case of — t 
is studied in the same way. 
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We let t = tq + e t, where t € [— to(a),to(a)] is a real-valued parameter. The 
number Ao is a double eigenvalue of %q{tq), and the perturbation e£o( T o + £t) is 
regular, see (|2 .2[) and (|2.3f) . and we can apply the standard regular perturbation 
theory, see e.g. [HI Ch. VH, §3]. The first of the estimates (|3.1I) shows that for 
sufficiently small £o > the operator H e (To + et) has two eigenvalues (counting mul- 
tiplicity) that converge to Ao as e — > +0. Denote these eigenvalues by E±(e,r) and 
the associated eigenfunctions by <fi±(x, r, e). Since the perturbation by eCofo + e t) 
is regular and the operators H e (to +e t) are self-adjoint, the eigenvalues E± (e, T+ei) 
and the associated eigenfunctions 0±(x, r + et, s) are holomorph w.r.t. e (the latter 
are holomorphic in Vt / | per (D)-norm), and the leading terms of their Taylor series 
near e = are of the form: 

E±(e, t + et) = A + eK±(t) + C(e 2 ), 

<p±(x,t,s) = V±(x,t) +e$±{x,t) +0(e 2 ), (3.3) 

:= cf* jtP + c+^ k , q , := c[¥jj, + C^^k.g, 

where cf — cf (t) are some constants that do not vanish simultaneously. Since the 
eigenfunctions <f>± can be chosen orthonormal in ^(D), the same is true for ^J±. 
This gives 

c+^ + c+^ = 0, \cf\ 2 + \cf\ 2 = l, (3.4) 

and the error estimates in (|3.3[) are uniform for t G [— to(a),to(a)] for any fixed a. 

To determine the coefficients K± it is sufficient to employ the regular perturba- 
tion theory. Namely, substitute the formulas (|3.3[) into the eigenvalue equation 

H e (t +et)<j>± = E±<j)± 

and equate the coefficients at the first power of s, then one arrives at the two 
equations 

(Hofo) - A )$± = 2<Z„(to)* ± + £(r )*± + K ± * ± . 

These equations are solvable iff the functions at the right-hand side are orthogonal 
to and # fc)9 in L 2 (D): 

(2« n (76)* ± + £(r )*± + K±^±, *j, P ) L2(D) = 0, 

(2tI n (7D)*± + £(r )*± + K ± *±, ^k, q ) L2(D) = 0. 

Using the orthonormality of \E" j. p and ^ k .q and Eqs. (|2.5[) and (|3.4[) we arrive at 

<^(ro)*j,p,*fc, 9 ) i2(n) = (ln(To)^k, q ,^ 3 , P ) L2{D) = 0, 

(Z„(r )% p , *j, p ) i2(D) = — — - t , (/„(ro)*fc,?, *fe,g) i2(D) = — y ~ T ° 



Thus, the solvability conditions (13.51) can be rewritten as 

(B(t) - K±(t))c±(t) = 0, B(t) = 2tB 1 - B (t ), 



,4\ f 2 -^ + r 



Therefore, K± are the eigenvalues of B(t) , and c± (t) are the associated eigenvectors. 
By an explicit analysis, 



K±(t) = h(r )t + k 2 (T ) ± yj (h(T Q )t + fc 4 (T )) 2 + |&? 2 (t )| 2 , (3.6) 

with fcj given by p7T2|) . By ([2~Tu]) one has #+(*) > for a11 * e K - Using the 

explicit expressions again we obtain 

2*1 (7b) = ^( T o) - -Qf-( T v)' 2fc 3(T ) - ^^(^o) - -gf-( T o)- 



8 



DENIS BORISOV AND KONSTANTIN PANKRASHKIN 



In accordance with (|2.8[) . the derivatives in the latter formulas have opposite signs, 
and this shows that 

|fci(r )| < |fc 3 (r )| ^ 0, (3.7) 
and implies the statement (IA1[) . Together with (|3.6[) it yields 

lim K±(t) = ±oo. (3.8) 

1 1 1 — >oo 

By elementary consideration one can find the minimum of K+ and the maximum 
of 

mm K + (t) = K+(t+), max K-{t) = K-(t-), 
* fe i( T o)| b i2( T o)| h(T Q ) 

*± = ± , = 7 — ) K±\P±) = P±(tq), (3.9) 

|Mi>)lv*§fo)-*?fo) 3 

where /3± are given by (|2.1ip . Let us choose the parameter a in (|3.2I) large enough 
and employ the asymptotics (13.31) . then 

min E + (e,r) = A + eK + (t + ) + 0(e 2 + er]), 

rG[0,7r/T] 

3.10 

max E-(e, r) = A + eKJtJ) + 0(e 2 + eri), 

Te[0.7r/T] 

and max E-(e.t) < min E-(e.t) for sufficiently small e due to (13.91) . 

re[0,7r/T] re[0,7r/T] 

Similar arguments are valid for r G [— 7r/T, 0], and one arrives at the analogues 
of (|3.3p and (|3.9p near — r . To obtain the analogues of the expressions (|2.13p . (|3.6p . 
(|3.7p . (13. 8p . (13. 9p and (|3.10p one should just replace (p,q,T ) by (— p, — g, — r ) in 
(|2~T2")) and (j2~TBl) . which gives exactly (|2~H)) and (j2~TBl) . This shows that He there 
exists a gap (a;(e), a r (e)) with a// r (e) given by (I2.18p . 

The extrema of K± are non-degenerate, \K±(t) — K±(t±)\ ^ C\t — t±\ 2 with a 
t-independent constant C, and 

\E±(e, t + et) - A - eisT±(t±)| > Ce|t - t±| 2 + C(e 2 + erf). 

This shows that the maximum points of E_ (e, •) and the minimum points of -E+ (e, ■) 
on t € [0,7r/T] are separated from tq + et + and tq + et_ by a distance of at most 
C(e 3 / 2 + erj 1 / 2 ), which proves the estimates (|2.19p and (|2.20p and completes the 
proof of the statement (|A2|) . 

To prove (|A3[) we remark first that = ^j.-p. If all the coefficients of the 
operator 7-L E are real-valued, then it commutes with the complex conjugation. In 
particular, C^{— T)^j- P = £o( r )^'i,p, an d /3±(tq) — /3±(— r ). Now the validity of 
(|2~T7|) follows trivially from ([AT]) . 

In the case To = the condition (|2.17p is reduced to the strict inequality /3_ (0) < 
/3+(0), which is true due to (lAip . Theorem 12. II is proved. 

4. Examples 

Let us discuss in greater detail several specific situations to which one can apply 
Theorem 12. II We will focus mainly on checking the conditions (|2.10p and (|2.I7D . 

4.1. Potential. Adding a real- valued periodic potential may be viewed as one of 
the simplest examples. One has simply C e — A§. In accordance with the conclusion 
(|A3P it is sufficient to check the condition (|2.10p . which loses the dependence on r 
and looks very simple, 
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This inequality is satisfied by a large class of potentials, and we can conclude that 
a generic potential perturbation opens the described gap in the spectrum as far as 
the cross-section and the period obey the relations discussed in Remark 12.21 



4.2. Magnetic field. Another natural perturbation is the action of a weak mag- 
netic field. The perturbed operator has the form W e = (iV + eA) 2 , where 
A = (Ai, . . . ,A n ) is magnetic vector potential, and the perturbing operator is 
given by 

*'' 9 -A^ + e\A\\ 



n , 
i=l ^ 



dxi dxj 



C (r)=iJ2[A i ^ + ^A i )-2rA, 



This results in 

n , 

rJ-r\ = i I A, 

' dxi dxi 

For t ^ we calculate the entries of the matrix Bq(t) using the integration by 
parts 



i=i 



L 2 (D) \ UJ "i> I L 2 (U 



2r(A n % p ,* fc , g ) L2(n) 
. n-l 



e ^^ At{x) (^-(x')Mx')- ^(x')Mx')) *» 
(^^ + r) Je 2 -^*~A n (x)^(x>)Mx')dx, 



6?i(r)=-|(^+r) jA n (x)^(x')dx, 

□ 

= - | + r) | A„(^(*') dx. 

□ 

Similar calculations for r < give the relations 

CM = -ffi i,m = 1,2, fci(-r) = -fci(r), i = 2,4, (4.2) 

and 

/3±(-t) = -^(r). (4.3) 
Taking into account (|Aip one can see that the condition (|2.17p becomes equivalent 
to /3+(t ) > ^_(-tq) and /3_(-t ) > /3+(r ). By (|P|) this reduces to 

±/3±(r )>0, (4.4) 

which is an equivalent compact form for (I2.17[) . 

Since p and q must have opposite signs (see Remark I2.2|) . one may assume 
without loss of generality p and g < 0, then the function ^(to) becomes 
positive. Substituting (|2.11l) . (|2.12|l and (|4~Tj) into one arrives at 

l^2(T-o)|^/fe|(To)-A:2(ro) > |fcifaM<r ) - fc 3 (ro)fc 2 (ro)|, 

(4.5) 



\bUro)\^kl(r )-kf(r )> 



2 \ 2 



/2?rp V f 27T< l 
an I — + T o I - ct22 I — + to 



where an := — / A n (x)^(x') dx and a 22 := ^ / ^(x)^^') dx. (4.6) 
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These inequalities together with 

b° 2 (T Q ) ± (4.7) 

give the final set of sufficient conditions guaranteeing the existence of a gap for T-L e . 

Let us now give an example of a specific magnetic field satisfying (14.5[) and (|4.7p . 
We restrict our attention to the functions A n (x) obeys 

T 



J A n (x)dx n — for all x' G ui, 



then (|4.6p shows that an = a 2 2 — 0, and the inequality (|4.5p is satisfied once the 
condition (|4. T|) is valid. Assume additionally A\(x) = A 2 (x) = . . . = A n -\{x) = 0, 
then the expression (|4.f [) for b\ 2 can be considerably simplified, and Eq. (|4.7[) 
becomes equivalent to 

T 

^(p + g) + ^ y j dXnAn{x)e ^i Xn ^ Q 

w 

For To = and t$ = tt/T the coefficient before the integral vanishes, see Remark l2.2l 
and our constructions do not allow to identify the gap opening (if a gap opens, its 
length is of order o{e)). On the other hand, for to G (0,it/T) the coefficient is 
non-zero, and all the required conditions are satisfied if one takes 

A n (x) = i}i(x')ip 2 {x')ipi{x') cos 27F ^ ^ x n + (p 2 (x), 

where ip\ G C^uJ) is an arbitrary positive function and tp 2 G C (O) is an arbitrary 
T-periodic w.r.t. x n function satisfying the condition 

T T 

I ifi2(x)dx n = / e T Xn ip 2 (x) dx n = for all x G u>. 



For instance, one take simply ipi = 1 and cp2 = 0. Thus we can conclude that the 
magnetic field can open a gap corresponding to the values of the quasi-momentum 
different from and ±7r/T. 

4.3. Deformation of boundary. Let us study a geometric perturbation consist- 
ing a periodic deformation of the boundary. We consider the two-dimensional strip 

■= {(2/1,2/2) : eh-(y 2 ) < 2/1 < l + eh + (y 2 )}, y = (2/1,2/2), 

where h± are smooth T-periodic functions. Denote by % £ the positive Dirichlet 
Laplacian in 51 e . Let us show that such operators are covered by the aforementioned 
construction. Let us denote h := h + — For small e the map $ : R 2 — > R 2 , 



, , (eh-{x 2 ) + (l + eh{x 2 ))x 1 

\XuX2j-y X2 



is the diffeomorphism, and f2 E = $ £ (£!o)- Let J$ £ be the Jacobi matrix of <£> e , then 
the map 

U e : L 2 (Q S ) ->■ L 2 (fi ), UJ = Vdet J$ £ /°$ £ , (4.8) 
is unitary, and the operator r H. e :— U e l-L e U* corresponds to the sesquilinear form 



&(«,«)=/ V( ^ ^ .(J^.J^vf = )detJ$ e (s) do:, 

Jn v Vdet J^e^r ^ ^/dct J$ e (x)^ 

(4.9) 
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i.e., 

1 



detJ$ £ V Vdet J$ e 

In our case 



, fl + eh(x 2 ) e(h'_{x 2 ) + x 1 h{x 2 ))\ 

■\J,%, X 2 ) — I q ^ I , 



which shows that "H e has real- valued coefficients and has the required representation 
H E = Ho + s£e ■ By direct calculation we get 



+ {h'_(x 2 )+ Xl h'(x 2 )){^ 



dx 2 dx 2 J dx\ dx\ 

du(x) dv{x) du(x) dv{x) 



dx. 



dx 2 dx 2 dx\ 

Let us find a class of functions h± for which the matrix entry (|2.10j) is non-zero. 
The transversal orthonormal eigenfunctions ipj are given by the explicit expressions 
ipj(x) = y/2 sin(7rja;), j = 1,2, Employing the notation 

h m = / e 2 * imx >' T h{x 2 )dx 2 , 



by direct calculations we obtain 

I := <£ (To)% P ,* fc , g ) L2(D) = (Co*^",*^")^ = - h + T +h , 
where 

□ 

h = 2 J h(x 2 )e 2 ^P-^/ T ^(xM Xl )dx = 2(wj)%- q 6 jk , 
□ 

^3 = i/ (/i'_(.T 2 ) + x 1 / l '( a;2 ))e 27ri ( p -^ 2 / T 
□ 

• (Jt + - jT-JV'j (ai^kfai) - (to + —)^ 3 (xi)ipk(xi)j dx. 

We observe that for j, k S N one has 

1 [0, j = fc, 

da* = (1 + { _ 1)j+k+1) _2jk . + ^ 
v j ~ k 

1 f -, j = k, 

x 1 ip j (x 1 )ip' k (x 1 )dx 1 =l 2 2 jk 

|(- 1 ) J+fe+1 -7 A p> 

Let us consider first the symmetric intersection To = 0. This implies immediately 
Q = ~P, 3 = k = 1, 



and 

1 /167r 2 p 2 „ ., . , 
I = -f- + 2n z h 



T V T 2 



i>2p- 



Thus, the condition ft,2p 7^ is sufficient to satisfy (|2.10[) . 



12 



DENIS BORISOV AND KONSTANTIN PANKRASHKIN 



For To = tt/T we have j = k = 1 and q = —p — 1 and 

1 = _ 1(433^ + *-))^, 

so the condition f|2.10[) holds true for fi2 P +i ^ 0. 

Let us consider the non-symmetric case To £ (0, tt/T). Without loss of generality 
we let j = 1 and k = 2. It implies immediately that I\ — I 2 — 0. On the other 
hand, by using (|2.21[) we get 

I 3 = -|(2t + 2^L±*>) jT (2hL<&) + h'(x 2 ))e^-^ T dz 2 

= 3 j? ( y + 2r J (ft- + Mp-g 

= 2n 2 (h77h + ) p _ q . 

Thus, the condition (j2~TU|) holds true as (ft_ + h+) p ^ 0. 

It is an interesting fact that for To = and To = n/T the sufficient conditions 
for a gap opening are formulated in terms of the Fourier coefficients of the function 
h := h + — h-, while for To £ (0,n/T) the same role is played by the function 
h + + h—. In other words, for To = and To = ir/T the gap opening is controlled, 
at the first order, by the strip width, while for t £ (0, tt/T) the same role is played 
by the sum of the side variations. 



4.4. Three-dimensional rod with a periodic twisting. Let 9 : K — > K be a 

smooth T-periodic function. For e > consider a diffeomorpisim $ e : M 3 — > R 3 , 

''cos (e9(x 3 ))xi — sin (e9(x 3 ))x 2 ' s 
® £ (x 1 ,x 2 ,x 3 ) = ( sin (e9(x 3 ))x 1 + cos (e9(x 3 ))x 2 

x 3 

and denote by ui we denote a two-dimensional connected domain with a piece- 
smooth Lipshitz boundary. Let Hq := us x M and il e :— $ £ (f2o)- For e = we 
just get a straight cylinder with a constant cross-section w, while for e 7^ the 
cross-section is rotated around the axis Oy 3 by the angle s9(y 3 ) w.r.t. the initial 
position in each plane y 3 = const. 

Let us denote by T-L e the Dirichlet Laplacian in f2 e . By analogy with (I4.9[) one can 
show that this geometric perturbation can be reduced to an additive perturbation. 
Since these constructions were already discussed in various contexts, we employ the 
ready expression for the perturbing operator obtained in [501 Eq. (15)]. Denote 
H s := U £ H e U* with the unitary U e given by (f¥75]t . then, for any u, v £ Wf (Q), 



(H s u, v} L2(n) = / V«(x) • G e {x)S7v{x) dx, (4.10) 



n 



where 



l + (ex 2 9'(x 3 )y -^xiX2(e , (x 3 )Y ex 2 9'(x 3 ) 
G e {x) = I -e 2 Xl x 2 (9'(x 3 )) 2 1 + (e Xl 9'{x 3 )) 2 -ex x 9'{x 3 ) 



ex 2 9'(x 3 ) —ex\Q'{x 3 ) 
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Thus, Ti. e = —A + sC £ , where the perturbation C s satisfies all the required condi- 
tions, and its leading term Cq is given by 



/ / On dv On dv \ 

(C u,v) L2{n) =J e'(x 3 )x 2 [—(x) — (x) + ——)(x) 



... . / du . . dv , , du , , dv , ,\ , 
-^)x 1 (_(x)_( a; ) + _(x)^(«))d af . 

We note that iH e commutes with the complex conjugation, and it remains to find 
the conditions guaranteeing the validity of (|2 . 10[) . There holds 

/ := <£ (ro)* J ,p,* fc , g ) L2(D) = (C^^e^ ^^e^)^ 

= i / e'(x 3 )x 2 e 2 ^-^/ T 
1 Jn 

. / 2nq \ dipj (2tt P \ dip k 

-i I — +to I ^^(^i, £2)% (£1,^2) +i I — + t I ^ii^i,^) g^-(xi,a; 

-0'(s 3 )^ie 27ri(p - 9)a;3/T 
. / 2irq \ dipj {2np \ dijj k 

= a(& — c), 

where 



d.r 



i| 9'( X3 )e^-^/ T d X3 = 2niP T2 ^ / «(.,,)<.--"■ ' <lr : ; 



& = i(^+ro) j /x 2 (^^ + ^ fc ) da* da* 
- i(-jr + to) _/ X! + ^ ) d* a dx 2 , 

c=(^±^ + 2r ) j^x 2 ^ k -x 1 ^ k )d Xl dx 2 . 

One observes that the coefficient a does not depend on the cross-section and can be 
chosen non-zero by an appropriate choice of the twisting function 9. The Green's 
formula gives the identities 

X2 (Vj i^r^ + Tr^^k) da;i dx 2 = [ x 2 -^—(i()jip k ) dxi dx 2 — [ x 2 t()ji) k dx 2 , 



~dx~ + 'dx'^ k ) dXl dX2 = j X2 'Q^r^ ) o^ k ) Axi ^ X2 = y x 2 ipjipkdx 2 

UJ 

- 1- - — Vfc d^i d a;2 = / siTr 

aa; 2 ox 2 / J da; 2 



Xl ^ j ~dx~ + ^x~^' k ) dxidx2 = f xi-^-(ipjipk)dxidx 2 = ~ f xiipjipk dxi. 



uj uj du 

Since the both functions t/^- and "0fc vanish at the boundary du>, one has 

xiipjipf. dxi — / x 2 ipj?pkdx 2 = 

J C?HJ dbJ 

which gives 6 = 0. 

Let us study the remaining coefficient c. In the case of a symmetric intersection 
of the unperturbed band functions, i.e., for To = or r = 7r/T, we have j = k and 
p q, which implies by (|2.21[) that c = and 7 = 0. Therefore, the assumptions 
of Theorem 12.11 are not satisfied, and the first order perturbation theory does not 
allow us to identify a gap opening for quasi-momenta close to the center and to the 
end points of the Brillouin zone. 
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For To ^ (0,7r/T) we have j ^ k and p ^ —q, and by (|2.2ip we conclude that 
p + q + 2to =/= 0. Let us show that the integral entering the expression for c does not 
vanish at least for some specific domains lu. Without loss of generality we let j = 1 
and k = 2, and take as an example lu = (0, n) x (0,7r/a) with a > 1. Denoting 
C = we have 

Ai = 1 + a 2 , ipi(xi, X2) = C sin(xi) sin(a!X2), 

A2 = 4 + a 2 , ip2(xi,X2) = C sin(2iEi) sin(aa;2), 

and 

-^-(x\, X2) = C cos(xi) sin(ax2), -^-(xi,X2) — aC sin(xi) cosfaa^)- 
0x1 0x2 

It yields 

dif}\ C 2 

x 2 -Q^-ip2 — — x 2 (sin(3xi) + sin(xi)) • (l — cos(2ax 2 )), 

dipi otC 2 , . . . .v . , . 

xi — — 1P2 = %i ( cos(xi) — cos(3xi)) ■ sm(2ax2 ), 

0x2 4 

and 

/ x 2 -^ — da;i dx 2 = — / (sin(3a;i) + sin(a;i)) dxi • / x 2 (l — cos(2aa;2)) da; 2 
J <JX\ 4 J Jo 

a 8 7T 2 4 
~ ^32? = 3a' 

7r Jr/a 

J xi— — ip2 dxi dx2 = ^ / si ( cos(xi) — cos(3a;i)) dxi J sin(2aa;2) dx2 = 0. 

cjj 

Finally, by employing (|2.21j) we obtain 

■ / d{x3)e 2 " i{p - q)x3/T dx 3 
./o 



/ 27r(p + g) , A 4 27r(p-g) 
i = —ac = — — h 2r n 



T 'J 3a T 2 

= -% f 9(x 3 )e 2 ^~^ T dx 3 , 
aT J 

and / ^ 0, if the Fourier coefficient in the previous expression is non-zero. 
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